function map2draw(lat,lon,raster,z)
% map2draw(lat,lon,image,z)
% -------- 自定义函数 --------
% 带经纬度绘制raster，使用pcolor
% lat：纬度向量
% lon：经度向量
% image: 是2维空间矩阵，现支持第三维是时间的三维矩阵，上北下南左西右东，
%   需要行数和lat大小一致，列数和lon大小一致。
% z：    可选参数，double，第三维的下标。
%
% 使用pcolor达到和image相同的绘制效果 
% pcolor绘制的时候是以左下为坐标起点的，使得绘制结果和image上下颠倒，用flipud
% 进行矩阵的倒置。
% 选用pcolor的原因是他不绘制nan数据，而imagesc会讲nan数据也给定颜色。

% 作者：sid(cnqdcyq@163.com)
% log：
%   2021/11/19 ver01 create

arguments
    lat
    lon
    raster {mustMatchCoor(lat,lon,raster)}
    z = 1
    
end

pcolor(lon,flip(lat),flipud(raster(:,:,z)));

shading flat
axis equal        
end

function mustMatchCoor(lat,lon,raster)
if any(size(raster)~=[numel(lat),numel(lon)])
    error('经纬度向量的大小和矩阵的大小不匹配！')
end
end